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ABSTRACT:  Isentropic  data  for  real  hydrogen  are  calculated  in 

a  range  of  pressures  up  to  30,000  atmospheres  and  temperatures 
up  to  3500°K.  The  effects  of  ionization,  excitation,  and 
dissociation  are  considered  negligible  because  of  the  rela¬ 
tively  low  temperatures  and  high  densities  involved.  The 
effects  of  the  intermolecular  force  are  accounted  for  by 
fitting  a  virial  coefficient  equation  to  low  temperature  high 
density  experimental  data.  These  calculated  data  in  turn  are 
fitted  to  an  empirical  entropic  equation  of  state.  This  par¬ 
ticular  equation  of  state  is  convenient  for  describing  many 
thermodynamic  processes. 
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CALCULATED  THERMODYNAMIC  PROPERTIES  OF  REAL  HYDROGEN  OP  TO 
30,000  ATMOSPHERES  AND  3500°K 

Modern  two-stage  gun  launchers  use  hydrogen  as  the  propellant 
gas.  The  prediction  of  the  performance  of  such  launchers 
requires  knowledge  of  the  hydrogen  gas  thermodynamic  data. 

This  paper  presents  calculated  thermodynamic  data  for  hydrogen. 
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INTRODUCTION 

Hydrogen  gas  is  often  used  as  the  propellant  in  light  gas 
guns.  In  the  so-called  "two  stage  guns"  the  hydrogen  propel¬ 
lant  may  be  compressed  to  densities  as  high  as  2,000  amagats*. 
The  temperature  is  increased  relatively  little,  going  from 
room  temperature  to  as  high  as  3500°K.  Under  these  conditions 
of  high  density,  large  intermoiecular  repulsive  forces  exist. 
The  hydrogen  gas  no  monger  can  be  described  by  the  ideal  gas 
equation. 

Unfortunately,  there  exists  no  reliable  experimental  data 
on  the  state  properties  of  hydrogen  at  high  pressures  and 
temperatures.  To  approximate  the  behavior  of  the  hydrogen  gas 
at  the  high  densities,  the  Abel-Noble  or  covoiume  equation  of 
state  is  sometimes  used.  This  equation  is 

p(v-b)  =  RT 


With  Cv  taken  to  be  a  constant,  the  isentrope  becomes 

p(v-b)^  =  K 


However,  the  Abel-Noble  equation  becomes  a  poor  approximation 
at  the  high  densities. 

Therefore,  a  calculation  of  hydrogen  gas  data  has  here 
been  made  which  attempts  to  better  account  for  the  inter- 
molecular  forces  due  to  high  density.  The  method  involves 
first  extending  the  calculation  of  Woolley,  Scott,  and 
Brickwedde,  reference  1,  for  the  thermal  properties  of  real 
hydrogen,  beyond  the  realm  of  experimental  data  as  reported  in 
reference  1.  Secondly,  the  resulting  theoretical  data  are 
fitted  to  the  semi-empirical  entropic  equation  of  Seigel, 
reference  2: 


p(8-2)/G 


(v-f) 


=  K 


where  B,  f,  and  K  are  constants  for  any  given  entropy.** 


*0ne  amagat  unit  is  the  value  of  density  at  1  atmosphere  and 
0°C. 

**This  equation  has  been  fitted  successfully  to  argon,  nitrogen, 

and  hydrogen  data  at  temperatures  below  430®K  and  pressures  up  «— } 

to  4,000  atm. 


t 
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It  is  seen  that  this  equation  is  equivalent  to  the  Abel- 
Noble  equation  for  a  given  entropy  if  3/(8 -2)  is  replaced  by 
and  if  f  is  replaced  by  b;  these  constants,  however,  will 
change  for  each  different  entropy.  The  semi-empirical  equa¬ 
tion  may  be  fitted  to  real  gas  data  with  much  greater  accuracy 
than  the  Abel-Noble  equation  because  of  the  fact  that  the 
parameters  B,  f,  and  K  may  vary  with  entropy.  In  addition, 
the  equation  is  convenient  to  use  to  describe  gas  dynamic 
processes  such  as  occur  in  the  two- stage  gun. 


A  BRIEF  DESCRIPTION  OF  REFERENCE  1 


In  reference  1  the  thermal  properties  of  hydrogen  are 
found  for  temperatures  up  to  5000°K,  but  at  low  densities 
where  molecular  interaction  is  negligible.  These  properties 
are  then  extended  to  somewhat  higher  densities  by  combining 
the  low-density  data  with  empirical  real-gas  data  at  inter¬ 
mediate  densities.  This  is  done  by  means  of  standard  thermo¬ 
dynamic  integrals  taken  along  isotherms. 

The  empirical  real-gas  data  of  reference  I  were  obtained 
in  the  following  manner.  The  empirical  equation 


Pv  -r{BP+CPZ) 
RT 


(1) 


was  used  to  account  for  the  intermolecular  forces.  The 
coeficients  B  and  C  were  assumed  to  be  functions  of  the 
temperature  only.  (When  the  above  exponential  is  expanded, 

B  is  found  to  be  the  second  virial  coefficient.)  The  effects 
of  excitation,  ionization,  and  dissociation  are  assumed 
negligible  because  of  the  relatively  high  density  and  low 
temperature.*  The  constants  B  and  C  were  fitted  to  experi¬ 
mental  hydrogen  data  in  the  realm  of  temperatures  from  273°K 
to  672°K  and  densities  up  to  500  amagat.  The  best  overall 
fit  by  a  weighted  method  of  least  squares  was  found  to  be  in 
reference  1: 

B  =  0.0055478T"1//i+-0.036877T"3//ij-.022004T”5//i| 

C  =  0. 004788t“3//2-0. 04053T~2 


*The  ionization  and  excitation  energies  are  very  high  (15.4  e.v 
and  13.5  e. v.,  respectively)  and,  hence,  the  effects  are  negligi 
ble;  the  dissociation  energy  (4.5  e.v.)  is  negligible  because 
of  the  high  densities  which  exist. 
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EXTENDING  TKE  CALCULATION  OF  REFERENCE  1 

In  figures  1  and  2  the  coefficients  B  and  C  are  plotted 
over  a  temperature  range  of  273°K  to  3500°K«  Reference  1 
carried  the  calculation  only  to  600°K,  for  this  is  the  extent 
of  the  realm  of  reliable  experimental  data.  Here  we  have 
extended  the  curves  for  these  coefficients  to  3500°K.  Values 
of  p,  v,  and  T  obtained  by  using  these  curves  with  equation  (1) 
are  listed  in  Table  1.  Values  of  S  calculated  by  the  method 
of  reference  1  are  listed  also. 

As  indicated  above,  the  effects  of  dissociation,  ioniza¬ 
tion,  and  excitation  of  the  hydrogen  are  assumed  negligible 
in  the  calculation.  This  is  probably  a  good  assumption, 
because  of  the  relatively  low  temperatures  and  high  densities 
considered.  The  Seigel  semi-empirical  equation 


p(B-2)/B  (v_f,  .  K 


(2) 


was  here  fitted  to  the  new  high-pressure  data  which  were 
calculated  above.  The  fitting  was  done  as  outlined  in 
reference  2  by  making  a  three-point  fit  of  equation  (2)  to 
the  data  along  a  given  isentrope,  thus  obtaining  the  empirical 
constants,  B ,  f ,  and  K.  Once  this  procedure  was  repeated  for 
a  number  of  isentropes,  the  constants  B ,  f ,  and  K  were  obtained 
as  functions  of  entropy  S.  Thus,  equation  (2)  was  fitted  to 
the  data  of  Table  1  at  temperatures  from  273°K  to  as  high  as 
3500°K  and  entropies  of  20,  22,  24,  26,  28,  30,  32,  34,  36, 

and  38  .  Figure  3  is  a  plot  of  the  P-V  data  of 

Table  1.  In  nearly  all  regions  of  this  plot,  agreement  between 
the  data  and  the  semi-empirical  fit  is  so  close  that  a  compari¬ 
son  cannot  be  shown  in  this  figure.  The  constants  b ,  f ,  and  K 
as  obtained  from  this  particular  three-point  fit  are  listed 
in  Table  2.  These  same  constants  are  plotted  against  entropy 
in  figures  4  through  6.  Other  thermodynamic  variables,  such 
as  internal  energy  U,  enthalpy  H,  or  sound  speed  a  can  be 
derived  from  equation  (2) .  The  following  expressions  (see 
reference  2  or  3)  result: 


U 


B-2 


+g 


H 


\  Kp2/fe  +f p+g 


B  v 

a  "  IPS  p  v=T 


r 
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where  g  is  a  function  of  S  only.  The  quantity  g  can  be  obtained 
in  a  manner  similar  to  that  of  obtaining  B,  f,  and  K.  Equation 
(2)  is  now  an  entropic  or  caloric  equation  of  state  from  which 
all  other  thermodynamic  variables  can  be  derived  in  standard 
fashion.  This  procedure  constitutes  a  rigorous  method  for 
describing  the  thermodynamics  of  a  particular  type  of  fluid. 

In  the  case  reported  here,  the  empirical  constants  of 
equation  (2)  have  been  fitted  to  data  in  a  region  lacking 
experimental  verification.  Though  this  means  of  obtaining 
real  hydrogen  data  does  not  have  the  desired  experimental 
foundation,  its  virtue  is  in  the  fact  that  it  continues  the 
trend  of  the  most  accurately  known  hydrogen  properties.  In 
all  likelihood  its  use  yields  a  better  approximation  than 
previous  techniques  used  to  predict  the  behavior  of  this  gas. 
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TABLE  I 


Properties  of  Hydrogen  from  Extended  Calculation  of 
Woolley,  Scott,  and  Brickwedde  of  NBS  (ref.  1) 


P0“14. 7psi. 
v0«178.365ft3/lb. 


u 

cal 

% 

V/v 

T 

mSIe^Tt 

Po 

vo 

(°K) 

20 

1066.42 

1.741  x  10”3 

300 

3000 

1.154  x  10"3 

395.570 

3130.29 

1.136  x  10~3 

400 

6000 

9.074  x  10~4 

473.285 

30,000 

5.617  x  10~4 

694.486 

22 

397.050 

3.464  x  10“3 

300 

1000 

2.135  x  10"3 

388.616 

1110.28 

2.029  x  10~3 

400 

3000 

1.304  x  10"3 

524.577 

6G00 

9.958  x  10“4 

632.060 

30,000 

5.856  x  10”4 

963.957 

24 

149.060 

8.037  x  10"3 

300 

408.842 

4.297  x  10“3 

400 

1000 

2.593  x  10~3 

514.004 

3000 

1.506  x  10"3 

696.100 

6000 

1.117  x  10“3 

840.715 

30,000 

6.212  x  10~4 

1289.86 

26 

55.4093 

2.047  x  10”2 

300 

-2 

100 

1.370  x  10 

355.417 

151.355 

1.038  x  10“2 

400 

1000 

3.202  x  10~3 

681.047 

3000 

1.776  x  10~3 

921.428 

6000 

1.277  x  10”3 

1109.88 

30,000 

6.677  x  10“4 

1682.07 
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TABLE 

S  P/ 

cal  /Po 

mole^K 

28  20.4089 

55.7425 
100 
1000 
3000 
6000 
30,000 

30  7.48250 

10 

20.4426 
100 
1000 
3000 
6000 
30 , 000 

32  2.73810 

7.48178 
10 
100 
1000 
3000 
6000 
30,000 

34  1.00124 

2.73602 
10 
100 
1000 
3000 
6000 
30,000 


(CONT.) 

V/vG  (3k) 


5.443  x  10““ 

300 

2.696  x  10“2 

400 

1.801  x  10“2 

472.335 

4.004  x  10“3 

901.583 

2.124  x  10~3 

1211.63 

1.480  x  10"3 

1449.26 

7.253  x  10“4 

2153.86 

1.473  x  10"1 

300 

1.200  x  lO"1 

326.090 

7.230  x  10“2 

400 

2.371  x  10"2 

627.815 

*■> 

5.038  x  10“° 

1186.74 

2.560  x  10~3 

1575.86 

1.731  x  10"3 

1868.39 

7.951  x  10"4 

2722.17 

4.015  x  10"1 

300 

1.963  x  10”1 

400 

1.597  x  10"1 

434.442 

3.123  x  10“2 

833.362 

6.336  x  10"3 

1545.34 

3.097  x  10”3 

2023.77 

2.036  x  10“3 

2379.57 

S. 788  x  10”4 

3406.82 

1.097 

300 

5.356  x  10"1 

400 

2.123  x  10”1 

578.167 

4.099  x  10”2 

1100.39 

7.934  x  10“3 

1986. 51 

3.748  x  10“3 

2568.45 

2.405  x  10”3 

2998.87 

9.785  x  10"4 

4229.22 
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TABLE  I  (CONT.) 


cal 

% 

V/ 

/  IT 

T 

5T^5K 

*  o 

V  O 

(°K) 

36 

0.36602 

2.999 

300 

1 

1.464 

399.974 

10 

2.819  x  10_1 

768.423 

100 

5.333  x  10“2 

1438.10 

1000 

9.878  x  10"3 

2522.81 

3000 

4.537  x  10“3  . 

3226.76 

6000 

2.850  x  10“3 

3745.12 

38 

0.13379 

8.204 

300 

0.36562 

4.003 

400 

1 

1.949 

532. 505 

• 

10 

3.728  x  10"1 

1016.94 

100 

6.857  x  10“2 

1855.11 

1000 

1.223  x  10~2 

3170.76 

3000 

5.486  x  10“3 

4017.95 

NOLTR  6 5-20 S 


£ 


TABLE  II 


Empirical 

Obtained 

Constants 
from  Three- 

of  Equation 
-Point  Fits 

(2)  vs. 
to  Curves 

Entropy  S  as 
of  Figure  3 

S 

B 

f 

K 

moTeOTT- 

vo 

P 

o  O 

20 

3.990 

2.81 

X 

10"4 

0.04735 

22 

4.  519 

3.15 

X 

1C'4 

0.08691 

24 

5.286 

3.44 

X 

10~4 

0.1677 

26 

5. 900 

3.65 

X 

10“4 

0.2798 

28 

6.338 

3.71 

X 

10""4 

0.4211 

30 

6.646 

3.65 

X 

10‘4 

0.5950 

32 

6.866 

3.41 

X 

10~4 

0.8166 

34 

7.025 

2C  60 

X 

io"4 

1.100 

36 

7.168 

4.00 

X 

10-5 

1.464 

38 

7.293 

-3.53 

X 

10“4 

1.949 

K 

1 


* 
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